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Abstract 

We study the Brodsky-Lepage evolution equation for the nucleon and 
construct an eigenfunction basis by including contributions of up to poly- 
nomial order 9. By exployting the permutation symmetry Pi 3 of these 
eigenfunctions, a basis of symmetrized Appell polynomials can be con- 
structed in which the diagonalization of the evolution kernel is consider- 
ably simplified. The anomalous dimensions are calculated and found to 
follow a power-law behavior. As an application, we consider the Brodsky- 
Huang-Lepage ansatz. An algorithm is developed to properly incorporate 
such higher order contributions in a systematic way. 

I. GENERAL FRAMEWORK 

The momentum scale dependence of the nucleon distribution amplitude is given 
by the Brodsky-Lepage evolution equation^ 



X1X2X3 
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(3 Jo 



[dy]V[xi,tfi]^(yi,0, (!) 



with £ = lnln-^ — ; Cb = = f an d Cf = %^ 1 = § the bosonic and fermionic 

Casimir-operators of SU(N) co i or , respectively; (3 — 11 — | Np = 9 being the Gell-Mann 
and Low function. The asymptotic solution of this equation is &as = I2OX1X2X3. To 
leading order in a s , the interaction kernel between quark pairs {i,j} is given in Ref. [I]. 
Factorization of Eq. (JTJ) leads to 



n xiX2X 3 <S>(xi,£,) = / [dy]V[xi,yi]$(yi,£) and xix 2 x 3 
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with 7 = f^f+^^f an d the integration measure is defined by J [dx] = Jq dxi Jq Xl d%2 

Jo 1 dx 3 5(1 — xi — X2 — £3). [In the following we use the operator Jo[dy] V[xi, yt] i— > V, 
which commutes with the permutation operator P13 [x\ £3).] The evolution be- 
havior of the eigenfunctions comes from the sensitivity on the transverse momentum 
integration^ arising from the gluon exchange kernels. It can be expressed in the form 

§(xi,Q 2 ) = ln~ 7F [jS — ) The representation of the evolution kernel is con- 
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veniently described in terms of Appell polynomials, ^F mn (xi,xs), which constitute an 
orthogonal basis with weights &As(xi) /120. It was shown in Ref. ^ that in this basis, 
V becomes block-diagonal with respect to different polynomial orders m + n. Because 
of [A3, V] = 0, it is usefulSi to define a symmetrized basis of Appell polynomials 

^Fmn(xi,x 3 ) = (1/2) (f mri (ii,i 3 ) ±J r nm (xi,x 3 )) for (m > n / m < n). (3) 

The particular importance of this basis lies in the fact that V is block-diagonal within 
a particular order for different symmetry classes (S n = ±1) with respect to Pi3- As a 
result, it is possible to analytically diagonalize V up to order 7 (cf. Tab. |) [up to order 
(0(n) =) M = 9 this is done in Refs. 3,4]. The eigenvalues of any order 0(n) and for 
a specific symmetry class S n form a "multiplet"-like set (see Fig. [l]), which follows an 
approximate power-law: 7„ = 0.37 0(n) ' 565 that differs from that of scalar (S=0) and 
vector (S=l) mesonsi'S. This observation is in contrast to the claims of Ref. [5[ The 
eigenfunctions of the evolution equation have properties of a commutative group 
00 

$ k (x l )®n(x l ) = J2 F LH^) with \0(k)-0{n)\ < 0(1) < 0(k)+0(n) (4) 
1=0 

and obey the above triangle relation. The structure coeflicientsi^ Fj. n of the group are 
calculated by F l kn = Ni ^[dx] x^x^il— x\— x s ) ^k(xi)^ n (xi)^i(xi), with the particularly 
important case F® k = j^. 

II. APPLICATIONS 

The evolution effect of the nucleon distribution amplitude is important for the 
calculation of various form factors at intermediate values of the momentum transfer 
Q 2 « 10 -h30 GeV 2 /c 2 (see Ref. UH-H). In this range, as{Q 2 ) tends to diverge whereas 
the Q 2 evolution of the distribution amplitude, given by 

B n (Q 2 )$ n (xi) , (5) 

n=0 / 

significantly reduces the value of the form factors by more than ps 30%. In the above 
representation of $tv, contributions of eigenfunctions up to polynomial order 3 have 
been studied&i. By incorporating on the rhs of eq. (||) (cf. Ref. 0) the factor f(xi, Xj) = 

e v 1=1 x i 2 J ; eigenfunctions of higher orders can be taken into account having 
recourse to a Brodsky-Huang-Lepage type of ansatz. This extended ansatz for Ai = 
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0.03 and A2 = 3 has been studied in Ref. |5|. In order to determine its evolution behavior, 
one has to project^ $at on the eigenfunctions $ n . For this purpose, / can be expanded 
in the basis of eigenfunctions of the nucleon evolution equation 

/fa, Xj) = ®k(xi), whh Ck (X j )=N k I [dx}x 1 x 2 x 3 f(x i ,X j )^ k (x i ). (6) 

k Jo 

This expansion and the properties of products of eigenfunctions given by the structure 
coefficients F^ k lead to 

$ N {xi) = ®As(xi) (^2 Bi Hx^ with B { = (^£B n c k F l n )j . (7) 

In the last equation the evolution of the distribution amplitude is fully determined by 
the scale-dependence of the expansion coefficients Bi(Q 2 ) = -E>z(/i 2 ) ( "^^j ) \ which 
are a complicated mixture of the original expansion coefficients B n and the projection 
coefficients c k of the extended ansatz. In contrast to the assumptionll 

^ n ( Xi ,q 2 ) = $ A5 (x,) ( n f^(p 2 ) (Spy) 7 " f{ - Xh Xj) > (8) 

the evolution of the distribution amplitude is not found to be determined by the scale- 
dependence B n (Q 2 ) = B n (/j, 2 ) ( ~^) ) 7 • The deviation of the approximation! from 
the correct evolution behavior can be expressed by the Q-dependent ratio 

TABLE I. Orthogonal eigenfunctions $ n (sci, x 2 , X3) = J2ik a ki x i x 3 °f the nucleon evo- 
lution equation (represented by the coefficient matrix a kl with a ki = S n a™ k ; a\\ = for all 



n). The normalization is given by J [dx] x\x 2 x 3 Q k (xi)& n (xi) = (N n ) 1 5 kn . 
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This deviation at intermediate Q 2 is found to be of order 0.9 for 0(7) < 3 and 
increases exponentially for higher orders. Since in the ansatz of Ref. [5] 70 > 7„ > 79, 
one can approximate 



FIG. 1. The eigenvalues of the evolu- FIG. 2. Deviation measure M.i{Q 2 = 30 
tion equation (for S n = 1) vs. the corre- GeV 2 /c 2 ) for the two extreme cases 7 = 70 and 
sponding order 0(n) = M (solid line). 7 = 79. 

Using the power- law behavior for ji, Fig. ^| shows the exponetially increasing 
deviation of the approximation of Ref. |5| compared to the Q 2 -scaling behavior based 
on the renormalization group equation. This comparison shows that it is important 
to project the nucleon distribution amplitude $at, as in the case of mesonsl, on the 
eigenf unctions of the evolution equation. For this purpose, powerful tools have been 
developed and the basis of eigenfunctions $ n has been extended up to polynomial order 
9. Using analytical and numerical algorithms developed in Refs. §,[|the calculation of 
higher order eigenfunctions up to any desired precision is shown to be possible. 
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